Ground-state properties of the two-band model for halogen-bridged metal complexes 
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Based on a symmetry argument, we systematically reveal Hartree-Fock broken-symmetry so- 
lutions of the one-dimensional two-band extended Peierls-Hubbard model. Performing numerical 
investigations as well, the possibility of novel density-wave states appearing is argued. 
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I. INTRODUCTION 



One-dimensional two-band models with competing 
electron-electron (el-el) and electron-phonon (el-ph) in- 
teractions have been attracting considerable interest. 
Such models cover quasi-one-dimensional materials such 
as halogen-bridged transition-metal (MX) linear-chain 
complexes and charge-transfer (CT) salts with mixed- 
stacks of alternating donor and acceptor molecules. 

One of the most interesting consequences of intrin- 
sic multiband effects and the competition between el- 
el and el-ph interactions may be the variety of ground 
states Mixed-stack CT compounds such as 

TTF-chloranil and TTeCiTTF-TCNQ may be described 
by a half-filled two-band model with site-off-diagonal 
el-ph coupling. They are dimerized into bond-order- 
wave (BOW) states at low temperatures, which is rec- 
ognized as a charge-transfer-induced spin-Peierls tran- 
sition. Ground states of MX compounds, which may 
be described by a 3/4-filled two-band model with site- 
diagonal el-ph coupling, strongly depend on the con- 
stituent metal and halogen ions. The conventional Pt 
compounds have the ground state with a charge density 
wave (CDW) on the metal sites induced by a dimeriza- 
tion of the halogen-ion sublattice, whereas recently syn- 
thesized Ni compounds P,|lO|] show a regular-chain struc- 
ture with a spin density wave (SDW) on the metal sites. 
Although the interchain hydrogen-bond networks are not 
negligible here, yet the competition between the CDW 
and SDW states may more or less be recognized as a 
crossover between the Peicrls and Mott insulators. Thus 
one-dimensional el-ph systems show rich phase diagrams 
and may possess further ground states unrevealed. 



II. MODEL HAMILTONIAN 



We study the one-dimensional two-band extended 
Peierls-Hubbard model: 
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Here we assume the Hamiltonian to describe MX chains. 



Now the notation of the Hamiltonian (2.1) is recognized 



as follows: 



,t 



and 



,t 



are the creation operators 



of an electron with spin s in the metal d z 2 and halogen 
p z orbitals at the unit cell n, respectively, u n the chain- 
direction displacement from the uniform lattice spacing 
of the halogen atom at the unit cell n, and L the number 
of unit cells; £m and £x are the energies of the metal d z 2 
and halogen p z orbitals on isolated atoms, respectively, t 
the transfer energy of hopping between these levels, Um, 
Ux, and V the el-el interactions on the metal atoms, on 
the halogen atoms, and between the neighboring metal 
and halogen atoms, respectively; a, /?, and K denote 
the site-off-diagonal and site-diagonal el-ph coupling con- 
stants and the elastic constant, respectively. We note 
that the model is rather general and may essentially be 
applied to various materials. For example, replacing the 
d z 2 and p z orbitals by pit orbitals of donor and acceptor 
molecules, respectively, we can immediately switch our 
argument to organic mixed-stac k C T compounds. We 
compactly rewrite Hamiltonian ( |2.l| ) as 
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where 
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are specified as 



(1 : k + q\t\l : k) = s M 5 q0 - -^ u *q sin (|) > 
(2 : k + q\t\2 : k) = £ x S q0 , 
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with the renormalized on-site affinities £m = £m — 
and e x = ex - U x /2. 

The symmetry group of the system is generally repre- 
sented as G = P x S x T, where P = Li A C2 is the 
space group of a linear chain with the one-dimensional 
translation group Li whose basis vector is the unit-cell 
translation, S the group of spin-rotation, and T the 
group of time reversal. Group actions on the creation 
operators are readily found in Ref. Let G denote 

the irreducible representations of G over the real num- 
ber field, where their representation space is spanned 
by the Hermitian operators {a\. k s a-j:k' ,s'}- There is a 
one-to-one correspondence fl2| between G and broken- 
symmetry phases with a definite ordering vector. Any 
representation G is obtained as a Kronecker product of 
the irreducible representations of P, S, and T over the 
real number field: G = P ® S ® T . P is characterized 
by an ordering vector q in the Brillouin zone and an ir- 
reducible representation of its little group P(q), and is 
therefore labeled qP(q). The relevant representations 
of S are S°(u(e,9)) = l and ^(ufe.fl)) = 0(u(e,6)), 
whereas ones of T are T°(t) — 1 T 1 (t) = — 1, where 
u(e, 9) = cr° cos(6>/2) — (<r-e) sin(0/2) represents the spin 
rotation of angle 9 around an axis e and 0(u(e, 9)) is the 
3x3 orthogonal matrix satisfying u(e,9)a x u'(e,9) — 
E^ x , y JO(u(e,9))] x ^ (A = x, y, z), with the 2 x 
2 unit matrix a and the Pauli-matrices vector er = 
(a x ,a y ,af). The representations P®5°(g)f , P®S 1 ®T 1 , 
Pg>S' <g)T 1 , and P<X>S ,1 <g)T correspond to charge-density- 
wave (CDW), spin-density- wave (SDW), charge-current- 
wave (CCW), and spin-current-wave (SCW) states, re- 
spectively. We leave out current-wave states in our 
argument, because here in one dimension all of them 
but one-way uniform-current states break the charge- or 
spin-conservation law. We consider two ordering vectors 
q = and q = tt, which are labeled r and X, respec- 
tively. Thus we treat the instabilities characterized as 



rP(r) <g> s° <g> f°, xP(x) ®s°® t°, rP(r) ® s 1 ® f\ 

and XP(X) ® S 1 <X> T 1 . Here P(r) and P(X) are either 
A (^-symmetric) or B (^-antisymmetric) representa- 
tion of C2 because P(P) = P(X) = C2 in the present 
system. 



III. BROKEN SYMMETRY SOLUTIONS 



In the HF approximation, the Hamiltonian (2.2) is re- 
placed by 
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(3.2) 
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where (• • -)hf denotes the quantum average in a HF 
eigenstate. 7Yhf can be decomposed into spatial- 
symmetry-definite components fl3| . Once a broken- 
symmetry Hamiltonian is given, its invariance group is 
defined. Keeping in mind that the density matrices char- 
acteristic of the given Hamiltonian should also be invari- 
ant for its invariance group, we can completely determine 
qualitative properties of the solution JTl] ] . 

We are interested in the local charge density on the 
site i at the nth unit cell, d i: „ = Yls( a t.n,s a i-n,s)BF, 
the local spin density on the site i at the nth unit cell, 

4.n = lY, s ,s'( a ln,s a i--n,s')iiF<J z ss ,, the complex bond or- 
der between the site i at the nth unit cell and the site 



2 



j at the rath unit cell, Pi-.n-j-.m = J2 s ( a L 
and the complex spin bond order between the site i at 
the nth unit cell and the site j at the mth unit cell, 
K.n-j-.m = \ Y< s ,s'(4-.n,s a 3--n,s>)wFO z ss ,. The halogen-atom 
displacements u n are self-consistently determined so as 
to minimize the HF energy _Ehf = (^)hf- All the solu- 
tions are explicitly written down elsewhere p3| , while we 
schematically show them in Fig. [j]. 



IV. NUMERICAL CALCULATION AND 
DISCUSSION 

We show in Fig. || phase diagrams at 3/4 band fill- 
ing, where q = tt states are predominantly stabilized. 
With the increase of £o, density waves on the halogen 
sites are generally reduced and single-band models come 
to be justified. The phase boundary between M-SDW 
and X-SDW is roughly given by Ux — Um — l-8eo, where 
£ = £m — £x- In order to realize X-SDW ground states, a 
small enough eo and a large enough Ux are necessary. We 
show in Fig. [| phase diagrams at half band filling, where 
q = states are relatively stabilized. The off-site el-ph 
coupling a stabilizes BOW, whereas the off-site Coulomb 
interaction V is unfavorable to that. 

Let us observe the doping dependence of the states. In 
Fig. ^ we plot macroscopic order parameters, Om-cdw = 

lEJ" 1 )"^:"' Om-SDW = 1 Z)n( _1 ) ns l:n' ^X-SDW = 
lEnf" 1 )"^) OboW = 2iEn (Pl:n;2:n-1 ~ Pl:n;2:n) , 
OsBOW = Yin (*l:n;2:n-l — *l:n;2:n); ^FM = 

and °FM = 7j£„ S 2:n. aS functions of band 

filling, where the parametrization is as follows: Uu/t = 
4.0, Ux/t = 2.0 for M-CDW, U D /t = 8.0, Ux/t = 3.0 for 
M-SDW, and Ur>/t = 3.0, Ux/t = 8.0 for X-SDW with 
common values e /t = 1.0, V/t = 1.0, a/(Kt)^ 2 = 0.1, 
and (3/(Kt)^ 2 = 1.0; U D /t = 6.0, Ux/t = 4.0, V/t = 1.0 
for BOW, U u /t = 10.0, Ux/t = 9.0, V/t = 2.0 for 
SBOW, and U D /t = 8.0, U x /t = 3.0, V/t = 1.0 for 
FM with common values e /t = 1.0, a/(Kt)^ 2 = 0.8, 
and (3/(Kt) 1/2 = 0.1. Now we explicitly find that X- 
phases are most stabilized at 3/4 band filling, while r 
ones arround half band filling. Though SBOW does not 
appear in the ground-state phase diagrams shown here, 
it seems to exist in the vicinity of half band filling with 
very strong on-site Coulomb repulsions. Even if the so- 
lutions are scarcely stabilized into a ground state, they 
can still be relevant in the ground-state correlations. We 
note, for example, that low-lying solitonic excitations in a 
SDW ground state induce SBOW domains around their 
centers [fL4"[ . We note in addition that 0^ M and Op M 
have opposite signs, which means spins on the metal and 
halogen sites are antiparallel to each other (Fig.|l|e). 

We finally stress the wide applicability of the present 
approach. Our approach never fails to reveal all the pos- 
sible broken-symmetry phases. We hope that the present 
argument will motivate further chemical, as well as the- 
oretical, explorations in low-dimensional el-ph systems. 
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FIG. 1. Schematic representation of possible density- wave 
states, where the variety of circles and segments qualita- 
tively represents the variation of local charge densities and 
bond orders, respectively, whereas the signs ± in circles and 
strips describe the alternation of local spin densities and 
spin bond orders, respectively. Circles shifted from their 
equidistant location qualitatively represent X-atom displace- 
ments. Identifying the present system, for example, with halo- 
gen-bridged metal complexes, each phase is characterized as 
follows: (a) Paramagnetism; (b) Electron-phonon bond order 
wave; (c) Metal charge density wave accompanied by an elec- 
tron-phonon bond order wave; (d) Halogen charge density 
wave accompanied by a purely electronic bond order wave; 
(e) Ferromagnetism with uniform spin bond orders; (f) Spin 
bond order wave; (g) Metal spin density wave accompanied 
by a spin bond order wave; (h) Halogen spin density wave 
accompanied by a spin bond order wave. 



FIG. 2. Ground-state phase diagrams at 3/4 band filling: 
(a) eo/to = 1.0 and (b) eo/to = 2.0 with common values, 
V/t = 1.0, a/(Kt) 1/2 = 0.1, and (3/(Kt) 1/2 = 1.0. 



FIG. 3. Ground-state phase diagrams at half band fill- 
ing: (a) V/t = 1.0 and (b) V/t = 2.0 with common values, 
eo/to = 1.0, a/(Kt)^ 2 = 0.5, and f3/(Kt)^ 2 = 0.1. 



FIG. 4. Order parameters as functions of band filling. 



3 



BOW 



M X M X M X M 
PM 



(a) 
(b) 
(c) 

(d) _ _^ 
^ X-CDW 

(e) ©se[+@+H+@+H+@ 

FM 





X-SDW 



Fig.l 
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